Paramagnetic Meissner effect in mesoscopic samples. 
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Using the non-linear Ginzburg-Landau (GL) theory, we 
study the magnetic response of different shaped samples in the 
field-cooled regime (FC). For high external magnetic fluxes, 
the conventional diamagnetic response under cooling down 
can be followed by the paramagnetic Meissner effect (PME). 
A second-order transition from a giant vortex state to a multi- 
vortex state, with the same vorticity, occurs at the second 
critical field which leads to the suppression of PME. 

PACS number(s): 74.24.Ha, 74.60.Ec, 73.20.Dx 

The Meissner effect is considered the most important 
characteristic property of superconductivity. When a su- 
perconductor is cooled down in the presence of an ex- 
ternal magnetic field the field is expelled and it behaves 
as a diamagnet. However, some samples show a param- 
agnetic response under cooling. The finding of PME (or 
Wohlleben effect) in high-T c superconductors initiated 
the appearance of several models interpreting PME as 
evidence of non-conventional superconductivity in these 
materials (e.g. see Ref. g]). However, numerous obser- 
vations of PME in conventional macroscopic l^-Q and 
mesoscopic || superconductors indicate the existence of 
another mechanism, which may be explained with GL 
theory. 

Based on results from axial-symmetric solutions of the 
GL equations by Fink and Presson ||, Cruz et al. || 
proposed that PME in their experiments on PbggTlgi 
cylinders is caused by a temperature variation of the su- 
perconducting density in a giant vortex state with fixed 
angular momentum J3|. In such a state, the supercon- 
ducting current, which shields the magnetic field in the 
vicinity of the sample boundary (essentially the Meiss- 
ner effect), changes its direction expelling magnetic field 
into the sample that can lead to the PME. Thereafter 
this idea, often called flux compression, was exploited in 
[||,|o|Q , but a quantitative analysis of PME within GL 
theory is still missing and several principal questions re- 
main to be answered: 1) is the vorticity of the giant vor- 
tex state fixed during cooling down as was assumed in 
Refs. 2) if yes, can this lead to the appearance of 

PME?, and, finally, 3) can the proposed mechanism ex- 
plain the PME in recent experiments with conventional 
macroscopic and mesoscopic || samples? In this 

Letter, we follow the GL approach and address these 
questions by studying the magnetic response of different- 
shaped samples in the FC regime. 



We consider a defect-free superconducting disk (and 
cylinder) immersed in an insulating media with a perpen- 
dicular (along cylinder axis) uniform magnetic field Hq. 
The behaviour of the superconductor is characterized by 
its radius R (and thickness d for the disk case) , magnetic 
field Ho, the coherence £(T) = £(0)(1 - T/Tc)- 1 / 2 and 
penetration A(T) = A(0)(1 - T/TJ- 1 / 2 lengths, where 
T c is the critical temperature. To reduce the number of 
independent variables we measure the distance in units 
of the sample radius, the vector potential A in ch/2eR, 
and the order parameter 'J in ^/—a/[3 with a, (3 being 
the GL coefficients |[2| . Then the GL equations become 

(-i% D -£f* = ^*{i (i) 

~A 3D A=^f(z)j 2D . (2) 

Here, the indices 2D, 3D refer to two-dimensional and 
three-dimensional operators; 

32D = ^ (**V 2Z5 4< - W 2 d**) - \M 2 A, (3) 

is the density of superconducting current in the plane 
(x,y), and the external magnetic field is directed along 
the z-axis. The boundary conditions to Eqs. (0-||) cor- 
respond to zero superconducting current at the sam- 
ple boundary and uniform magnetic field far from the 
sample. The order parameter and superconducting cur- 
rent are assumed to be independent of the z— coordinate 
which is valid for cylinders as well as for thin « 
disks gf-|lj|. Then f(z) = 1 and f(z) = d5{z) for the 
cylindrical and disk geometry, respectively. 

Superconducting disks and cylinders placed in a mag- 
netic field and cooled down transit from a normal to 
a superconducting state at the critical temperature 
which depends both on Hq and R. However, the unitless 
parameters Hq/ H^T*), R/^T-A, and the angular mo- 
mentum of the giant vortex superconducting state in the 
nucleation point T* depend only on the magnetic flux 
<I> = itR 2 Hq piercing through the sample Here, 
H C 2 = < f>o/27r^ 2 and <f>o = hc/2e are the second critical 
field and the flux quantum, respectively. With further 
cooling down, the magnetic response is characterized by 
only two independent variables Hq/H C 2 and 

Our numerical approach for solving Eqs. (0-||) is de- 
scribed in p3[-p^|. It turns out that an accurate simu- 
lation of a multi-vortex state is a hard task in the case 
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of large vorticity L. The latter correponds to the total 
angular momentum and the number of vortices in the 
giant-vortex and the multi-vortex state, respectively. To 
improve the accuracy we apply a non-uniform rectangu- 
lar space grid condensing in the vicinity of the sample 
boundary. However, due to the tremendous computa- 
tional expenses (e.g. the total number of grid points in 
the plane (x,y) was about 160000 for L = 30) we could 
only treat the vortex state with large L in the cylindri- 
cal case, when the vector potential is uniform in the z- 
direction. In the disk case, we restrict our considerations 
to the axial symmetric solutions \P = ?p(p)exp(iL(j)) (p, <p 
are the cylindrical coordinates), which are shown to be 
stable in the region Hq > H C 2- 

When the sample is cooled down below the critical 
temperature a giant vortex state appears with angu- 
lar momentum L, which is determined by the magnetic 
flux $ H^|Q. Starting from this state we mimick the 
FC regime by decreasing (increasing) slowly the value of 
H /H c2 ~ 1/(1 - T/T c ) (R 2 /e ~ (1 - T/T c )) such that 
the system evolves along a path with fixed external mag- 
netic flux (e.g. see Fig. [I]). Using the superconducting 
state found at the previous step as input, we find the 
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next steady-state solution to Eqs. (|I[ 
consider only stable solutions and neg 
tuations, which could lead to possible transitions between 
metastable states. This assumption is valid for nor- 
mal superconductors where the barriers separating the 
metastable states exceed by far the sample temperature, 
except near points in which the state becomes unstable 
[ p7[ , e.g. near the saddle points. 

When calculating the dipole magnetic moment, we 
can neglect non-linear effects in the vicinity of the nu- 
cleation point. The quantum angular momentum L in- 
creases almost proportional to the magnetic flux $ but 
remains always smaller than $/$o- The supervelocity 
v,f, = p~ l {L — $p 2 /&oR 2 ), which is oriented along the 
azimuthal direction, changes its sign at p+ = R^L^o/®. 
Therefore, both diamagnetic (p > p*) and paramagnetic 
(p < p^) currents exist in any giant vortex state. How- 
ever, the magnetic moment, which can be estimated in 
the lowest Landau level (LLL) approximation as D ~ 
J dpp v^ipLf with ipz being the lowest eigenfunction of 
the linearized first GL equation, turns out to be always 
diamagnetic both for disks and cylinders. As long as the 
superconducting density |^| 2 remains small, the mag- 
netic moment almost linearly increases in absolute value 
with decreasing Hq/H C 2 (inset (a) in Fig. [l]). With fur- 
ther cooling down, the LLL approximation breaks down. 
Due to non-linear effects (mainly from the second term in 
the RHS of Eq. ([!])) the order parameter increases more 
rapidly in the inner region which leads to an increase of 
the paramagnetic component (Figs. |l].||). 

The resulting magnetic moment crucially depends on 
the ratio between L and With increasing magnetic 

field, the switching L — > L + 1 of the angular momentum 
of the nucleated state occurs at a certain [[14]. For 



given L, the magnetic moment reaches its maximum and 
minimum value at and respectively. Due to 

angular momentum quantization the magnetic moment 
exhibits a strong oscillating behaviour as function of the 
magnetic field (see inset of Fig. |J) which agrees with 
Geim's observations ||. With decreasing temperature, 
the dipole magnetic moment becomes zero for certain ra- 
tio H/H C 2 ~ 1/(1 — T/T c ) and the original diamagnetic 
response becomes paramagnetic. This magnetic field is 
shown in Fig. ||for $ = and <I> = $£,. Note that: 1) 

for the cylinder geometry a larger GL parameter k = A/£ 
favours PME, 2) while an increase of the effective pene- 
tration length X 2 /d suppresses PME in disks, and 3) an 
increase of $ favours PME both in cylinders and disks. 
The reason is that the point p*, where the supervelocity 
changes its direction, shifts towards the sample bound- 
ary with increasing angular momentum which is related 
to the total magnetic flux $(L > 1) w $ (£ + £ 1/2 ) @- 

In cylinders, the magnetization is directly proportional 
to the magnetic moment. In disks: 1) the diamagnetic 
currents flowing near the sample boundary give larger 
contributions to the magnetization than inner paramag- 
netic currents [jllj and consequently, the paramagnetic 
contribution to the magnetization will be strongly sup- 
pressed (inset (b) in Fig. Q), 2) but the smaller trapped 
magnetic flux, in disks as compared to cylinders, de- 
creases the diamagnetic response. This is the reason wh 
thicker disks show onset of PME at larger Hq / H C 2 (Fig 
open symbols). 

When the second critical field H C 2 (T) becomes smaller 
than the applied magnetic field, the giant-vortex state 
transits to the multi-vortex state with the same vortic- 
ity (Figs. |^,^). This second-order transition is not fol- 
lowed by any jumps in the magnetization or the mag- 
netic moment. Just after the transition, all vortices are 
arranged in a ring (0:L). Note, that the magnetic moment 
of the state (0:L) continues to increase with decreasing 
H a /H c2 {T) but with a smaller slope (Fig. g). With fur- 
ther decreasing Hq/H C 2(T) a pair of vortices moves to 
the inner region and the state (2:L-2) appears for L = 20 
(see Fig. |^). This first-order transition is followed by 
a weak jump in the magnetic moment. The derivative 
dD/dT changes sign and further cooling down results 
in the disappearance of PME. The magnetic moment of 
the diamagnetic state with smaller angular momentum 
L = 19 and $ = is also affected by the transition to 
the multi-vortex state, which increases the diamagentic 
response (inset (b) in Fig. |J). As the temperature de- 
creases, vortices continue to move from the outer to the 
inner shell. Note, that the corresponding weak jumps in 
the magnetic moment are practically not visible on the 
scale used in Fig. [|. Although the state with L = 19 is 
energetically more favourable than the one with L = 20, 
no vortex exits the system which is in agreement with 
observations ||. Since the vorticity remains unchanged 
under cooling down, the magnetic response will be al- 
most reversible. The hysteresis caused by the first-order 
transitions between different multi-vortex states with the 
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same vorticity is weak (see inset (a) of Fig. [|). Starting 
from the point H / H C 2 = 0.9 and warming up the system 
to H/H C 2 = 1.1 we find the magnetic moment, which co- 
incides in the scale of Fig. ^ with that obtained by cooling 
down. 

Although the magnetization curves, shown in Fig. ^, 
agree qualitatively with those from experiments a 
number of issues remain unclear: 1) in experimental ob- 
servations of macroscopic samples a weak hysteresis is 
found in cooling down and subsequent warming up, and 
2) the observed maximum in the magnetic momentum 
are weaker than those from our calculations which may 
be due to the presence of vortex pinning centra. Note 
that within GL theory we found a weak PME which is 
caused by the competition of large diamagnetic and para- 
magnetic responses of the outer and inner part of the 
sample, respectively. Any mechanism which slightly in- 
fluences any of the two responses may strongly influence 
the total magnetic behavior. As an example, in macro- 
scopic disks PME disappears after mechanical abrading 
the top and bottom surfaces A number of sam- 

ples made from the same material as those demonstrating 
PME exhibited only diamagnetic behaviour |I| ||. This 
indicates the important role played by the sample struc- 
tural inhomogeneity. For mesoscopic disks, experiments 
pj show PME for rather small angular momenta which 
does not agree with our simulations for flat circular meso- 
scopic disks. 

To address some of these sample structural issues we 
consider the influence of the superconductor shape on 
the magnetic moment by varying radially the thickness 
of the disk. We limit ourselves to the case of a strong 
type-// superconductor (k S> 1) and solve, therefore, 
only Eq. ([!]). The dipole magnetic moment of different- 
shaped samples is shown in Fig. (||) for H = H c2 (T), 
where the scale of the thickness variation is apparent 
from the insets of Fig. ^. The magnetic moment becomes 
more negative (positive) in magnifying glass (crown) like 
samples as compared to flat disks. An increase of the 
local thickness near the sample boundary (see Fig. ||(c)) 
increases the PME. This strongly suggests a possible non 
flat geometry of the disks of Ref. || . 

In summary, the giant-vortex state remains stable un- 
der cooling down and transits to the multi-vortex state 
with the same vorticity at Hq w H C 2- The paramag- 
netic response is caused by a more rapid growth of the 
superconducting electron density in the inner region of 
the sample, due to non-linear effects, where paramag- 
netic currents flow. The appearance of the multi-vortex 
state supresses PME and the maximum of the paramag- 
netic response corresponds to Ho rj H C 2. We showed that 
within the GL theory a paramagnetic response is possi- 
ble for large magnetic fluxes, which is in agreement with 
experimental findings on disks with large radia, but does 
not agree with the experimental results on mesoscopic 
disks of Geim et al ||. 

After finishing this work, we came aware of a preprint 
by Palacios |ll| who used the LLL approximation to 



study the FC regime and failed to find any paramag- 
netic response. As shown above, one has to go beyond 
the LLL approximation in order to find PME which is 
caused by non-linear effects. 
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FIG. 1. The nucleation field as a function of the ra- 
dius. Thin curves show the path along which the system 
evolves during cooling down. Insets (a) and (b) show the 
corresponding dipole magnetic moment and magnetization 
(M = J dr(H z - H )/dr, wh< ire integration is performed 
over the sample volume), respectively. Here, solid and dot- 
ted curves correspond to cylinder (k = 10) and disk (k = 1, 
d — 0.1R), respectively. 

FIG. 2. The distribution of the superconducting electrons 
density (a) , the superconducting current (inset) and the mag- 
netic field (b) in the cylinder (solid and dashed curves) and 
the disk (inset) for different reduced magnetic fields H/H C 2 : 
1-1.6, 2-1.3, 3-1.0. 

FIG. 3. The magnetic field corresponding to the transition 
from the diamagnetic to paramagnetic response as a function 
of the total angular momentum for <E> = and $ = <3>l. 

Solid and open symbols correspond to cylinders and disks, 
respectively. Inset shows the magnetic moment of the cylin- 
drical sample at H — H C 2 and k = 10. 

FIG. 4. The dipole magnetic moment of cylindrical sam- 
ples for different external magnetic fluxes $ = and 
k = 10. The hysteretic behaviour is demonstrated in inset 
(a) for L = 20 and $ = Inset (b) shows the magnetic 

response for L — 19 and $l- 

FIG. 5. The dipole magnetic moment and contour plots of 
the superconducting density (log-scale is used) in the vicinity 
of the second critical field. 

FIG. 6. The dipole magnetic moment of different shaped 
circular samples at Hq/H C 2 = 1. 
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